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Abstract 

We discuss initial conditions for the recently proposed Imperfect Dark Matter (Mod¬ 
ified Dust). We show that they are adiabatic under fairly moderate assumptions about 
the cosmological evolution of the Universe at the relevant times. 


1 Introduction and Summary 


Planck data favors adiabatic initial conditions at the onset of the hot Big Bang [T]. This 
corresponds to the situation in the early Universe, when the number densities of different 
particle species, e.g., dark matter (DM) particles and photons, are universally distributed in 
space. That picture, comfortably accomodated in the ACDM cosmology supplemented by 
the short stage of an inflationary expansion, can be less transparent in more exotic setups. In 
particular, adiabaticity of initial conditions is obscured in the case, if DM has a non-particle 
origin. 

In the present paper, we continue to discuss Imperfect Dark Matter (IDM) scenario [H |3] 
(Modihed Dust in Ref. [2]). The action of IDM is given by 


S 



^ - 1 ) + 


( 1 ) 


(see Ref. [1] for generalizations). Here A and ip are two scalars, and 7((^) is some function 
of the held ip. The Lagrange multiplier A enforces the constraint g^'^df^ipd^ip = 1. This 
equation dehnes the geodesics followed by the collisionless particles, with the held ip being 
the velocity potential. The analogy with the pressureless perfect huid, dust, is complete in 
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the situation when the higher derivative (HD) term is absent, i.e., 7 ((^) = 0. In that case, 
the Lagrange multiplier A plays the role of the energy density of dust: it redshifts away 
as the inverse third power of the scale factor a in the expanding Universe. While the dust 
provides a good description of DM to the linear level, it should be modihed in the non-linear 
regime. The reason is that it develops caustic singularities, i.e., the physical quantities,—the 
velocity dispersion and the energy density,—blow up at the finite time mu- This drawback 
of the dust model was one of the motivations (not the unique one, though) to introduce the 
HD term [2j. Before we summarize some of the effects arising due to the non-zero function 
7 (</?), let us discuss, how the action ([T]) emerges in different physical frameworks. 

The mimetic dark matter scenario mi deals with the conformally transformed metric 
of the fornji]. 

Rather unexpectedly, variation of the general relativity action with respect to the auxilliary 
metric gap and the scalar gives more than just Einstein’s equations. The difference is 
about a new degree of freedom behaving as dust. This returns us to the action ([1]) up to 
the HD term [10], which is set by hands in this context [ 8 ]. Alternatively, one can trace 
back the origin of IDM to Lorentz violating theories of gravity. In particular, the action ([T]) 
arises in the infrared limit of the Horava-Lifshitz model with the projectability condition 
applied mi da da]. More generically, IDM is closely related to a version of Einstein- 
Aether theory [T3| with the aether held being the derivative of the scalar mds] (see also 
Refs. midHldSI for the most recent discussion on the topic). 

During the large part of IDM evolution, one assumes that the function 7 ((^) is constant. 
In this situation, the model possesses shift symmetry ^ tp + c. There is, consequently, the 
Noether charge density that redshifts away as the inverse third power of the scale factor. Up 
to the term suppressed by the 7 factor, the energy density of IDM is equal to the Noether 
charge density [3| (Section 2). Hence, with a good accuracy, the cosmological evolution of 
IDM mimics that of the dust particles B The degeneracy gets broken at the linear level, 
where the 7 -term leads to the constant sound speed ~ 7 [Sj. This sets a cutoff on the power 
spectrum at sufficiently small scales: beyond the sound horizon energy density perturbations 
do not grow. Therefore, they are suppressed compared to the predictions of cold dark 
matter (CDM) scenarios. In particular, setting]^ 7 ~ 10“®, one can suppress the growth of 
structures with the comoving wavelength < 100 kpc, alleviating the mismatch between the 
observed number of dwarf galaxies and the value predicted in the CDM framework [SU] Ell 

^The generalization of mimetic dark matter scenario to disformal metrics transformations has been con¬ 
structed in Ref. |9]. 

^This is an exact statement in the situation, when IDM is the only component of the Universe. See the 
discussion in Section 2 

•^Note that we use the convention SirMp^ = 1 throughout the paper, where Mpi is the Planck mass. An 
associated value of the parameter 7 reads in units of Ref. m 7 - 
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[ 22 ] , Alternatively, however, the small scale problems may be the consequence of the incorrect 
implementation of several baryonic processes, as indicated most recently in the analysis of 
Ref. [23]. In that case, one rather deals with the constraint [2] 

7 < 10 -^ . ( 2 ) 

For those small values of the parameter 7 , the linear evolution of IDM perturbations is 
analogous to that of CDM given that they start from the same initial conditions (see the 
discussion below). The important difference, however, may arise at the structure formation 
level, i.e., in the non-linear regime [2]. This line of discussion is far out of the scope of 
the present paper. Here we will be interested in the opposite situation, namely, when the 
relevant cosmological modes are in the deep super-horizon regime. 

If the shift symmetry is exact at all the times, IDM cannot be the main component of 
the invisible matter in the Universe. The reason is that the Noether charge density gets 
washed out during inflation with an exponential accuracy. In this situation, IDM constitutes 
only a tiny, 0 ( 7 ), fraction of the overall DM during the dust dominated stage [3]. To avoid 
this, one necessarily assumes breaking of the shift symmetry taking place at the early stages 
of the Universe, i.e., deeply in the radiation dominated (RD) era. In IDM model, this is 
realized by promoting the constant 7 to the function 7 (</?) 0 . We assume furthermore that 
the variation of the function is substantial only at very early times, and negligible 

otherwise. In that way, one can easily generate the amount of the Noether charge required 
to explain cosmological experiments, as we review in Section 2. Notably, perturbations of 
IDM produced by the same mechanism are adiabatic with a high accuracy. 

Before the Noether charge gets produced, IDM tracks the dominant matter of the Uni¬ 
verse, e.g., it has an equation of state of radiation in the Universe driven by the relativistic 
particles. At the linear level, the similarity with the dominant matter reveals in an exact 
adiabaticity of IDM perturbations (Section 3)@. That situation is quite analogous to what 
one has in the ACDM cosmology at very high temperatures, when all the particle species are 
in the thermodynamic equilibrium. More importantly, the adiabaticity of IDM perturbations 
holds later on, after the shift-symmetry breaking occurs, despite the presence of the non- 
adiabatic pressure. We prove this by the explicit computation of the curvature perturbation 
of IDM, i.e., Cidm- We show that the latter corresponds to the adiabatic initial conditions for 
IDM, at least under rather moderate assumptions: the variation of the gravitational poten¬ 
tial must be negligible relatively to the Hubble rate at the times, when the shift-symmetry 

"‘Another idea would be to couple the field y to the inflaton [7]. In the present paper, however, we consider 
only gravitational interactions between IDM and other fields. 

®The analogous observation has been made in the earlier work (2). Here we show that this is an exact 
statement. We also generalize it to any dominant matter in the Universe, while the case of radiation has 
been considered in Ref. [2]. 
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breaking takes place. That condition is obeyed with a high accnracy well before the niat- 
ter/radiation eqnality. In particular, this guarantees that IDM is indistinguishable from 
CDM at the level of cosmic microwave background measurement^. 

The outline of the paper is as follows. In Section 2, we briefly review the IDM scenario 
including the mechanism for producing the Noether charge density in that picture. In Sec¬ 
tion 3, we discuss the super-horizon evolution of IDM perturbations and show that they are 
adiabatic under rather general conditions. 


2 Generating dark matter 

We start with writing down the system of equations following from the action ([1]). The 
simplest one is the constraint 

= 1 , (3) 

obtained from the variation with respect to the Lagrange multiplier. Applying the covariant 
derivative to the constraint ([3]), one obtains the geodesic equation. In this regard, IDM is 
equivalent to the collection of dust particles moving in the gravitational held. This degen¬ 
eracy, however, gets broken at the level of the other equations. Variation of the action ([1]) 
with respect to the held ip yields [5] 


= :p^(apf , 


where 


Finally, the energy-momentum tensor corresponding to the action ([1]) is given by mn 


(4) 

(5) 


Ti: =i^X- 27 ^ 

- 7((p) . 


+ 7 ( 9 ?) ( 




In the present paper we are mostly interested in the background evolution of IDM as well 
as in the linear evolution of its super-horizon perturbations. That is, we will always neglect 
spatial derivatives of the helds. In this approximation, Eq. dl]) takes the forrrEI, 




( 6 ) 


®This way of setting initial conditions is to be compared to that of Ref. [2]. In the latter paper, IDM 
perturbations relied on the arbitrary constant of integration. That is, the adiabaticity of initial conditions 
in Ref. [2] was at the price of tuning an arbitrary constant to some particular value. Alternatively, any other 
choice would result in an admixture of an isocurvature mode. 

^We stick to the notion ’’Noether charge”, that can be misleading at the times, when the shift symmetry 
is broken. Hopefully, this is not going to confuse the reader. 
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where n = is the Noether charge density. We assume that during most of the stages of 
the evolution, the function is constant, i.e., the model is shift symmetric. In that case, 
the r.h.s. of Eq. (|4]) is zero and one has 


d{na?) 

dt 


0 


n = X — S'jH . 


(7) 


As it follows, the Noether charge density redshifts away as 1/a^, i.e., n = C/a^. Assuming 
that IDM obeys shift symmetry during the inflationary expansion of the Universe, the density 
n relaxes to zero with an exponential accuracy, i.e., n = 0 at the onset of the hot Big Bang. 
Otherwise, one would need to set an extremely large value of n in the beginning of inflation. 

In the situation with the zeroth Noether charge, IDM tracks the total matter hlling in 
the Universe [3]. Indeed, the physical energy density is related to the Noether charge density 
by 

37 

PiDM = Tq = n + —ptot , ( 8 ) 

while the pressure is given by 


V,dm = -ij’* = . 


(9) 


Here we made use of equations = \ptot and —2H — 3H^ = Vtot] the subscript HoP stands 
for the total matter including IDM. Neglecting the Noether charge density as well as the 
variation of the parameter 7 in Eqs. (jH]) and ([9]), one obtains wjdm = where 

Wtot = ^ 7 ' This proves the point made above. 

For the future convenience, we prefer to rephrase Eqs. ([ 8 ]) and ([9]) by splitting the con¬ 
tributions to the total energy density and pressure, which follow from IDM and the external 
matter helds. This gives 


2n 37 

Pi DM 7^ 7, 1 “ T, o dext 

2 — 37 2 — 37 


( 10 ) 


and 

Ddm ^ . (11) 

Here p^xt = Ptot - Pidm and Text = T’tot - T’/dm denote the total energy density and 
pressure of all the matter helds not including IDM, respectively. As it follows from Eqs. flTU]) 
and flTT]) . IDM behaves as a perfect tracker also with respect to the external matter. Namely, 
setting the Noether charge density and the derivative of the parameter 7 to zero, we have 
widm = uJext, where w^xt = In particular, this means that for a very small value of 

pext 

the parameter 7 , IDM may constitute only a tiny fraction of the invisible matter during the 
dust dominated stage. We are interested in the different opportunity of IDM being the main 
(the only) component of DM. 
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Since this point on and until the end of the Section, we assume that there is a period of 
shift-symmetry breaking taking place deeply in the RD stage. This allows to generate non¬ 
zero Noether charge, which, by the suitable choice of parameters, can be tuned to match the 
observed value [5]. With no much loss of generality, we consider an instantaneous transition 
from some initial value 71 to the value 72 > 71 - Namely, the function 7 ( 95 ) has the form 

7(‘^) = 7i + A7 ct((p - (p*) , (12) 

where (T((p — (p*) is the Heaviside function; (p* is the constant pointing the time, when the 
transition happens, and A 7 = 72 — 71 . 

Accordingly, the value of the Noether charge undergoes an instantaneous flip. Integrating 
Eq. ([2]) with the initial condition n = 0 set at t < t*, we obtain 

9 

na^ = 2^7 ' \t=u ^ ^ (^ 2 ) 

One can tune the time and the quantity A 7 , so that we get the correct energy density of DM 
today. As it follows, earlier the Noether charge is generated, less value of A 7 is required. 
In particular, for the choice A 7 ^ 7 10“®, this indeed happens at very early times. The 

corresponding redshifts and temperatures read .2 > 10^^ — 10^^ and T > 100 MeV — 1 GeV, 
respectively. Soon after the Noether charge is produced, it becomes the total contribution to 
the energy density of IDM. Indeed, the second term in Eq. ffTOj) mimicking the behaviour of 
the external matter of the Universe (radiation) redshifts away fast relatively to the Noether 
charge density. Since this point on, the cosmological evolution of IDM resembles that of the 
dust particles. This is an exact statement in the situation, when IDM is the only matter in 
the Universe. In that case, pidm = ptot and Vjdm = Vtot- From Eqs. (IH]) and ([9]), one then 
easily obtains pi dm oc n oc 1 /a^ and Vjdm = 0 [2], [3l ITS] . 

In the remainder of the paper, we show that perturbations of IDM generated by the same 
mechanism are adiabatic under fairly relaxed assumptions. 

3 Super-horizon evolution of IDM perturbations 

3.1 Generalities 

Before we dig into the details of the linear level analysis, let us specify the gauge choice. 
Generically, the metric reads to the linear order, 

= (1 + 2$)(dx°)^ -I- 2adiZdx^dx^ — a^{6ij + 2'^dij — didjE)dx^dx^ , 

where E, Z, $ and d' are the scalar potentials, and we omitted vector and tensor perturba¬ 
tions. See Refs. [2ll [25l [26] for the reviews and textbooks on the cosmological perturbation 
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theory. Calculations look particularly simple and transparent in the synchronous gauge, 
which we use in the bulk of the paper. In the Appendix A, we provide calculations in the 
Newton’s gauge in order to cross-check the results. 

To go from the Newton’s gauge, where Z = 0 and Z = 0, to the synchronous gauge, one 
makes the coordinate transformation, + where obey the condition 

doii + diio = 0. With this coordinate choice, the perturbation 5ip and the potential $ turn 
into zero, 

5ip ^ 5(p = Q $—)-|) = 0. 

The fact that perturbations of the held and the potential <h can be switched to zero 
simultaneously follows from the gauge-independent relation 50 = <h,—the constraint ([3]) 
linearized. In the synchronous gauge, the potential T is given by 

^ ^ ^ - H5ip . (14) 

The condition di^Q + do^i = 0 then guarantees that the (Oi)-component of the metric remains 
zero. At the same time, the potential E is generically non-zero in the synchronous gauge. 
This is, however, negligible in the super-horizon regime, which is the case of our primary 
interest. Indeed, dtdjE = —{di^j + dj^i)/a? oc didjSip/{a^H) —)■ 0. 

Hereafter, we prefer to omit the tilde over the transformed quantities. In the synchronous 
gauge, the Noether charge density perturbation is given by (see Eq. (]5|)), 

5n = 5\ — SyT — SyT . (15) 

In what follows, we will also need the expressions for IDM energy density and pressure 
perturbations. These are given by 


^PiDM = 6n + 9'yH'^ 


( 16 ) 


and 

5Ei£)m ~ — 9'yH^ — SyT — SyT , (f^^) 

respectively. 


3.2 Before shift-symmetry breaking 

We first discuss the case of the zeroth Noether charge. Let us show that perturbations of 
IDM are exactly adiabatic in that case. One obtains from Eq. flTbl) . 

SpiDM = . (18) 
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The analogous expression for the energy density perturbation of the external matter can be 
inferred from the 00-th component of Einstein’s equations, which reads in the synchronous 
gauge 

1 1 

(19) 


Pi dm 


37 


-5p, 


ext 


Combining this with Eq. flTSl) . one obtains 

iWDM = 2 _ 3 ^. 

We plug the latter into the definition of the curvature perturbation of IDM, 

SpiDM 


CiDM = + 

and then make use of the relation 

PiDM + Vidm = 


^{PlDM + 'Pidm) 

37 


( 20 ) 


2-37 


i^Pext T ext) ■ 


This expression follows from Eqs. (nsii and dj, where one should set the Noether charge 
density to zero and the parameter 7 to a constant. We obtain, 


ClDM — d/ -I- 


6p, 


ext 


3{p 

ext + 'Pext) 


— Cext ) 


( 21 ) 


where C,ext is the curvature perturbation of the external matter. Finally, we use the relation 
between the partial curvature perturbations corresponding to the external matter helds and 
the quantity (ext, 

(ext= 

i^IDM 

Here the index i stands for the particular matter field, i.e., photons, neutrinos, baryons etc., 
and the subscript i 7 ^ IDM means that the contribution of IDM has been omitted. Under 
an assumption that there is no admixture of the baryon or neutrino isocurvature modes, all 
the partial curvature perturbations Q^idm are equal between each other. Since an equality 
Y.i^iDMPi/Pe.xt = 1, one has = (ext- Comparing this with Eq. ([21]), we get 


(iDM — (i^IDM ■ 


( 22 ) 


That is, the IDM curvature perturbation equals to the curvature perturbations of the stan¬ 
dard matter fields. This means that IDM perturbations are exactly adiabatic in the shift- 
symmetric case/before the shift-symmetry breaking takes place. That situation is quite 
similar to what happens in the standard cosmology with the particle DM: at very early 
times all the particle species behave as the single fluid, and initial scalar perturbations are 
adiabatic by definition. 






3.3 After shift-symmetry breaking 

While IDM perturbations are exactly adiabatic before shift-symmetry breaking, it is not 
immediately clear that they remain so at later times. Indeed, starting from initial condi¬ 
tions (12^ . they can change due to the presence of the non-adiabatic pressure. 


= SVlDM - , (23) 

PiDM 

which is not manifestly zero. This is one distinction of the IDM fluid from the familiar fluids, 
i.e., radiation and dust. Consequently, the IDM curvature perturbation evolves behind the 
horizon accordingly to the equation |24] . 


ClDM 


H 

PiDM + VidM 


q^non—ad 
' IDM 5 


(24) 


while those of the standard matter fields remain constant. This results into the violation of 
the condition fl2^ leading to the appearance of an isocurvature mode and, consequently, to a 
potential conflict with the Planck data. The non-adiabatic pressure, however, is negligible in 
two regimes: at the times t <t^ and t ^ t*. The former follows from an exact adiabaticity 
of IDM perturbations at the early times implied by Eq. fl22lfl The latter is also clear, since 
IDM relaxes to the standard dust at sufficiently late times. Hence, the non-adiabatic pressure 
is relevant only in the intermediate regime t ~ , when the effects due to the shift-symmetry 

breaking may become strong enough. They are encoded in the appearing Noether charge 
density in Eqs. (IHI) and (|T6|) and explicitly in terms involving the derivative of the function 
7 in Eqs. (|9]) and (|T71) . These new terms source the non-adiabatic pressure. To summarize, 
an expected change induced in the curvature perturbation C/dm is measured in terms of the 
quantities calculated at the times t ~ t*. Shortly, we will confirm this observation by making 
an exact calculation. We will also see that the variation of the curvature perturbation Cidm 
is small, as it relies on the derivative of the potential T. 

To study the evolution of the curvature perturbation, Eq. fl2TD is not very convenient. 
For this purpose, it is simpler to exploit Eq. (|6]). Integrating the latter and using Eq. flT^ . 
one obtains at the times t > t^, 


(5n = 3n(T(t,)-T) + 2n-. (25) 


®Not referring to Eq. (l2^ . one can show that the non-adiabatic pressure of IDM at the times t is 
proportional to that of the external matter, i.e., oc where the subscript ’ext’ stands for the 

combination of photons, neutrinos, baryons etc. As perturbations of the external matter fields are assumed 
to be adiabatic, one has = 0. Consequently, P77 ia 7°^ = should be. 
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Using then Eqs. (E]), (jl]), (fT 6 |) . (fT^ and (1^ . we derive the expression for the curvature 

perturbation of IDM, 


3n(^(U) - ^) + 6 n • ||^ + 
3[?7, — ?>'yH] 


ClDM — fl/ + 


(26) 


Hereafter, 7 denotes the value of the parameter at the end of the transition 71 72 , i.e., 

7 = 72 - Note that at the times t < t*, when the Noether charge density equals to zero, 
Eq. fl26|) reduces to the expression for the curvature perturbation of the total matter (tot, as 
it should be. Now, we are interested in the different regime, when IDM mimics the behaviour 
of the dust particles, i.e., strong inequalities n ^ 'j\H\ ~ 'jH^ ~ '’jptot are obeyed. For the 
relevant values of the parameter 7 this still happens deeply in the RD stage. In that limit, 
we get for the curvature perturbation of IDM, 


(iDM — ^'(U) + 2 


H{Q ■ 


(27) 


The latter is constant as expected: the non-adiabatic pressure is negligible in the late-time 
regime. Now, let us show that the second term on the r.h.s. of Eq. fl27)) equals to zero, 
i.e., the derivative of the gravitational potential vanishes at f = f*. This we will do in the 
reasonable approximation, when all the external matter is in the state of radiatioij^. We will 
need the ij-component of Einstein’s equations, which reads in the synchronous gauge, 

T -|- = —— {ST’ext + 5'Pidm) ■ (28) 

For the relativistic external matter, one has 5Vext = \^Pext- Using this, we combine Eqs. flT^ 
and fl25]) to exclude the quantities describing the external matter, 

T -|- 477\i/ = -6piDM — -^5 Vidm ■ 

6 2 

We substitute expressions flTO]) and (ITT)) for the IDM energy density and pressure perturba¬ 
tions and rewrite the equation above as follows. 


d 

dt 




+ 477 1 



6n 


Integrating this out, we get 


1 - 


37 


= 


a^(i)5n{i)di + C 


(29) 


®In particular, this statement is exact for the times t* corresponding to the temperatures T > 100 GeV, 
when all the Standard Model degrees of freedom are relativistic. 
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Here ti < t and C are the arbitrary constants. Recall now that the Noether charge density 
is zero at the times t < t^. As we are interested in the behaviour of the potential tl' at t 
and since the quantity n is always hnite, the hrst term on the r.h.s. of Eq. fl29|) vanishes. 
Hence, the solution for the derivative of the gravitational potential reduces to 





(30) 


For the arbitrary value of the constant C, this solution is discontinuous, as it follows from the 
behaviour of the parameter 7 at f = f*. Note, however, that the r.h.s. of Eq. fl30l) translates 
into the decaying mode of the potential which is commonly dropped in cosmologj0. That 
is, we should set the constant C to zero. In this situation, the derivative of the potential 
T vanishes, i.e., = 0, and the expression for the IDM curvature perturbation flTTI) 

simplihes to (idm = d/(f*). 

Consequently, the potential T is a continuous constant function at the time t = t^. From 
Eq. fl25|) . we then conclude that 6n{t^,) = 0. The same is true for the IDM and external 
matter energy density perturbations, i.e., SpjoMit*) = 0 and 5pextii*) = 0. These follow 
from Eqs. fllbp and 0191) . respectively. In the situation, when all the degrees of freedom are 
relativistic and there is no admixture of the isocurvature mode in the particle sector, one 
has an equality 5ext = ^ph-, where 6ext = and 6ph = 7 ^; the subscript "ph" stands for 
the photons. Hence, = 0. Recall now the expression for the curvature perturbation 

of the photons, 

Cph d' + ■ 

The latter stays constant behind the horizon, i.e., Cph(t) = Cp/i(C)- Thus, (ph = d/(t*). To 
summarize, there is no admixture of the IDM isocurvature mode, i.e.. 


SlDM,ph = ^{(iDM — (ph) — 0 . (31) 

Perturbations of IDM remain adiabatic after the shift-symmetry breaking, provided only 
that we can neglect the non-relativistic degrees of freedom during the transition 71 —?■ 72 . 

The result fl3T]) can be also understood from a slightly different prospective. As we noted 
earlier, the IDM curvature perturbation starts from exact adiabatic initial conditions 
and may change only due to the non-zero adiabatic pressure fl23ll . The latter can be relevant 
only at the times f ~ t*. Next, we observe that both the energy density perturbation (the 
numerator in Eq. fl26|) I and the pressure perturbation flTT)) . rely only on the variation of the 
potential d/. Accordingly to Eq. ([23]), so does the non-adiabatic pressure and, consequently, 

^°Recall that the scale factor during the RD stage grows as a oc ^/i. Hence, the decaying mode of the 
potential d/ drops as 4'dec oc y. 
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the isocurvature perturbation. Hence, the resultant perturbations must be adiabatic, once 
the potential \1/ is constant at the relevant times. 

Finally, let us comment on the generality of the results obtained. First, we notice that 
the assumption of the instantaneous transition 71 —)■ 72 is not a strong one at all. Our results 
hold for the different choices of the function 7 (v 5 ) given that the gravitational potential \ 1 / 
remains constant during the phase of the shift-symmetry breaking. That condition is satisfied 
with a high accuracy provided that the transition occurs well within the RD stage. 

Second, instead of the HD term as in Eq. ([1]), one could consider another one, 

+ • (32) 

We provide the associated analysis in Appendix B. In particular, we show that all the three 
statements take place: i) in the shift-symmetric case, the equation of state of IDM is that 
of the total matter; 2) perturbations of IDM are exactly adiabatic in this situation; 3) an 
approximate adiabaticity holds after the short phase of shift symmetry breaking. 

To conclude, initial conditions for IDM are the same as in the CDM case. Hence, CDM 
and IDM result with the same predictions regarding the cosmological observations. An 
important difference, however, may arise at the galaxy scales. This is due to the fact that 
IDM possesses non-zero sound speed. Furthermore, IDM and CDM exhibit an apparently 
different behaviour in the non-linear regime [2] . We leave this and other interesting questions 
for the future. 
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Appendix A 

In this Appendix, we switch to the Newton’s gauge in order to cross-check the results obtained 
in the main body of the paper. In that case, the perturbation of the Noether charge density 
is given by 

Sn = SX + 6'yH^ + 3'yH^ — 37 T , 

where the shift-symmetry is assumed. The expression for the energy density perturbation of 
IDM reads in the Newton’s gauge, 

SpiDM = 5n — 9777 ^$ -|- 9'yH'^ . (33) 

Recall that the Noether charge density equals to zero before the shift-symmetry breaking 
takes place. Integrating Eq. ([6]) with the initial condition n = 0 set at the times t < one 
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obtains at the times t > U, 


6n = -2n ■ <l>(t*) + 3n ■ ■ VtotiQ + 2n ■ , (34) 

where n ■ a? = const is defined by Eq. ffT^ . Substituting Eq. fl54)) into Eq. fl55]) . and then 
the latter into fl2U]) . we obtain for the IDM curvature perturbation, 


ClDM — 4^ + 


1 


3[n — 37Ef] 


X 


3n(4'(t*) - 4^) + n • (t*) + ^7 ■ 


■VtotiU) , (35) 


where 5tot = Sptot/ptot- Here we made use of the 00th component of Einstein’s equations, 
which reads in the Newton’s gauge (in the super-horizon regime), 

3H^ - 3 / 7 ^$ = . 

The curvature perturbation fl35D corresponds to exactly adiabatic initial conditions for IDM 
before the Noether charge is produced, as expected. Well after the shift-symmetry gets 
broken, i.e., in the regime n ^ one has 


C/DM — 4/(f*) -|- + 


1 

3 H{t^) 


'Ptotit*) ■ 


This expression is equal to the curvature perturbation (|27)) calculated in the synchronous 
gauge, as it should be. In particular, neglecting the super-horizon variation of the potential 
T at the time f* and using dip ~ and H ^ P, we get 


ClDM ~ 4/(t*) -I- 


^ Ptotii*) 
4ptot (^*) 


~ Cph • 


This corresponds to the adiabatic initial conditions for IDM. 


Appendix B 


In this Appendix we consider another possible HD term given by Eq. fl52]) . The associated 
energy momentum tensor is given by [15] 

T^u = Xd^ipdyip + ^V^iVxipV^'Vup) - ^DV^ipVuP - . 

At this level, we assumed the shift-symmetry, as it is going to be enough for our purposes. 

First, it is straightforward to show that in an exactly shift-symmetric case, IDM still 
tracks the dominant matter of the Universe. Indeed, from the energy-momentum tensor one 
deduces for the energy density and the pressure. 


PiDM — n + —ptot 


(36) 
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and 


'PiDM — '^'Ptot ■ ( 37 ) 

Here n is the Noether charge density now given by the expression 

n = J^ = Xd\ - V" (yV.VV) 

As it follows, in the situation, when the Noether charge is zero, the equation of state of IDM 
is that of the total matter in the Universe. Not surprisingly thus, super-horizon modes of 
IDM are exactly of the adiabatic type. 

To show this explicitly, we set to zero the Noether charge density n as well as its super¬ 
horizon perturbation, 

5 n = 5 \ + 6777\i/ . 

The relation between the IDM energy density perturbation and the quantity 5 n is given by 

5piDM = + . ( 38 ) 


Hence, 

SpiDM = SyTTii/ . 

Substituting this into the dehnition (1^ . using Eqs. fl36D and (1^ . one obtains Cidm = Ctot- 
The latter implies an exact adiabaticity of IDM perturbations in the situation, when the 
Noether charge density equals to zero. 

To obtain the expression for the perturbation C/dm in the generic case, one exploits an 
equation, 

h . „) = . (39) 

Namely, we promote the constant 7 to the function of the held (p. We again assume with no 
loss of generality, the instantaneous transition of the parameter 7 as in Eq. (IT^ . Integrating 
Eq. (I39|) over the time, we get 

3 

n ■ t > . 

The analogous equation at the level of super-horizon huctuations coincides with Eq. (|25|) . 
Substituting Eq. (]25|) into Eq. (l38l) . we obtain for the IDM energy density perturbation 

SpiDM = 3n ■ (^'(U) - ^') n ■ 6tot{U) + ■ 6ptot ■ (40) 

We substitute this into Eq. fl2U]) . At sufficiently late times, i.e., when the strong inequality 
n ^ ^ptot is obeyed, we arrive at the expression flTT)) . Again neglecting the super-horizon 
variation of the potential d' and then following the same steps as in the bulk of the paper, 
one concludes with adibaticity of IDM perturbations. 
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